d(H2O) / d(t) = (-2)*k1*em^2 + (-1)*k4*em*H + k6*H*OH + (-2)*k9*em*O2m + k11*H*H2O2 + k15*OH*H2O2 + k17*OH*HO2 + (-1)*k19*HO2*O2m + (-2)*k20*O2m^2 + k21*H2O2 + k22*Hp*OHm + (-1)*k23*H2O + (-1)*k26*HO2m + k27*OHm*H2O2 + (-1)*k30*Om + k31*OH*OHm + (-1)*k34*O2m + k35*OHm*HO2 + (-1)*k38*em + k39*H*OHm + k40*OH*H2 + (-1)*k41*H + k42*OH*HO2m + PH2O
d(em) / d(t) = (-2)*k1*em^2 + (-1)*k4*em*H + (-1)*k5*em*OH + (-1)*k7*em*H2O2 + (-1)*k8*em*O2 + (-1)*k9*em*O2m + (-1)*k10*em*HO2 + (-1)*k36*em*Hp + k37*H + (-1)*k38*em + k39*H*OHm + Pem
d(Hp) / d(t) = (-1)*k22*Hp*OHm + k23*H2O + (-1)*k24*Hp*HO2m + k25*H2O2 + (-1)*k28*Hp*Om + k29*OH + (-1)*k32*Hp*O2m + k33*HO2 + (-1)*k36*em*Hp + k37*H + PHp
d(H) / d(t) = (-2)*k2*H^2 + (-1)*k4*em*H + (-1)*k6*H*OH + (-1)*k11*H*H2O2 + (-1)*k12*H*O2 + (-1)*k13*H*HO2 + (-1)*k14*H*O2m + k36*em*Hp + (-1)*k37*H + k38*em + (-1)*k39*H*OHm + k40*OH*H2 + (-1)*k41*H + k45*H2*Om + PH
d(OH) / d(t) = (-2)*k3*OH^2 + (-1)*k5*em*OH + (-1)*k6*H*OH + k7*em*H2O2 + k11*H*H2O2 + 2*k13*H*HO2 + (-1)*k15*OH*H2O2 + (-1)*k16*OH*O2m + (-1)*k17*OH*HO2 + k28*Hp*Om + (-1)*k29*OH + k30*Om + (-1)*k31*OH*OHm + (-1)*k40*OH*H2 + k41*H + (-1)*k42*OH*HO2m + POH
d(OHm) / d(t) = 2*k1*em^2 + k4*em*H + k5*em*OH + k7*em*H2O2 + 2*k9*em*O2m + k16*OH*O2m + k19*HO2*O2m + 2*k20*O2m^2 + (-1)*k22*Hp*OHm + k23*H2O + k26*HO2m + (-1)*k27*OHm*H2O2 + k30*Om + (-1)*k31*OH*OHm + k34*O2m + (-1)*k35*OHm*HO2 + k38*em + (-1)*k39*H*OHm + k43*H2O2*Om + k44*HO2m*Om + k45*H2*Om + POHm
d(H2) / d(t) = k1*em^2 + k2*H^2 + k4*em*H + (-1)*k40*OH*H2 + k41*H + (-1)*k45*H2*Om + PH2
d(H2O2) / d(t) = k3*OH^2 + (-1)*k7*em*H2O2 + k9*em*O2m + (-1)*k11*H*H2O2 + (-1)*k15*OH*H2O2 + k18*HO2^2 + k19*HO2*O2m + k20*O2m^2 + (-1)*k21*H2O2 + k24*Hp*HO2m + (-1)*k25*H2O2 + k26*HO2m + (-1)*k27*OHm*H2O2 + (-1)*k43*H2O2*Om + PH2O2
d(HO2) / d(t) = (-1)*k10*em*HO2 + k12*H*O2 + (-1)*k13*H*HO2 + k15*OH*H2O2 + (-1)*k17*OH*HO2 + (-2)*k18*HO2^2 + (-1)*k19*HO2*O2m + k32*Hp*O2m + (-1)*k33*HO2 + k34*O2m + (-1)*k35*OHm*HO2 + k43*H2O2*Om + PHO2
d(HO2m) / d(t) = k10*em*HO2 + k14*H*O2m + (-1)*k24*Hp*HO2m + k25*H2O2 + (-1)*k26*HO2m + k27*OHm*H2O2 + (-1)*k42*OH*HO2m + (-1)*k44*HO2m*Om + PHO2m
d(Om) / d(t) = (-1)*k28*Hp*Om + k29*OH + (-1)*k30*Om + k31*OH*OHm + (-1)*k43*H2O2*Om + (-1)*k44*HO2m*Om + (-1)*k45*H2*Om + k46*O3m + (-1)*k47*Om*O2 + POm
d(O2) / d(t) = (-1)*k8*em*O2 + (-1)*k12*H*O2 + k16*OH*O2m + k17*OH*HO2 + k18*HO2^2 + k19*HO2*O2m + k20*O2m^2 + 0.500*k21*H2O2 + k46*O3m + (-1)*k47*Om*O2 + PO2
d(O2m) / d(t) = k8*em*O2 + (-1)*k9*em*O2m + (-1)*k14*H*O2m + (-1)*k16*OH*O2m + (-1)*k19*HO2*O2m + (-2)*k20*O2m^2 + (-1)*k32*Hp*O2m + k33*HO2 + (-1)*k34*O2m + k35*OHm*HO2 + k42*OH*HO2m + k44*HO2m*Om + PO2m
d(O3m) / d(t) = (-1)*k46*O3m + k47*Om*O2 + PO3m

k1 = 5.950e+06
k2 = 1.040e+11
k3 = 9.860e+09
k4 = 4.960e+11
k5 = 3.720e+11
k6 = 6.340e+10
k7 = 2.850e+11
k8 = 2.180e+11
k9 = 1.610e+11
k10 = 1.610e+11
k11 = 2.140e+09
k12 = 6.060e+10
k13 = 2.130e+11
k14 = 2.130e+11
k15 = 4.230e+08
k16 = 8.970e+10
k17 = 3.200e+10
k18 = 4.100e+07
k19 = 4.810e+08
k20 = 3.500e+07
k21 = 3.780e-03
k22 = 1.130e+12
k23 = 6.520e-02
k24 = 5.690e+11
k25 = 2.520e+01
k26 = 1.760e+08
k27 = 1.360e+11
k28 = 5.690e+11
k29 = 2.520e+01
k30 = 1.760e+08
k31 = 1.360e+11
k32 = 5.690e+11
k33 = 1.550e+05
k34 = 2.870e+04
k35 = 1.360e+11
k36 = 7.160e+11
k37 = 1.650e+05
k38 = 2.010e+03
k39 = 8.010e+09
k40 = 7.800e+08
k41 = 2.100e+03
k42 = 8.180e+10
k43 = 8.180e+10
k44 = 8.760e+10
k45 = 1.550e+09
k46 = 1.990e+07
k47 = 3.260e+10

G_nH2O = 0.000e+00
G_nem = 1.290e+00
G_nHp = 1.290e+00
G_nH = 5.900e-01
G_nOH = 2.870e+00
G_nOHm = 0.000e+00
G_nH2 = 1.000e+00
G_nH2O2 = 4.000e-01
G_nHO2 = 3.000e-02
G_nHO2m = 0.000e+00
G_nOm = 0.000e+00
G_nO2 = 0.000e+00
G_nO2m = 0.000e+00
G_nO3m = 0.000e+00

G_gH2O = 0.000e+00
G_gem = 3.430e+00
G_gHp = 3.430e+00
G_gH = 1.560e+00
G_gOH = 5.740e+00
G_gOHm = 0.000e+00
G_gH2 = 6.400e-01
G_gH2O2 = 2.700e-01
G_gHO2 = 0.000e+00
G_gHO2m = 0.000e+00
G_gOm = 0.000e+00
G_gO2 = 0.000e+00
G_gO2m = 0.000e+00
G_gO3m = 0.000e+00

Pem = 1.0364E-7*rho*(D_g*G_gem+D_n*G_nem)
PHp = 1.0364E-7*rho*(D_g*G_gHp+D_n*G_nHp)
POH = 1.0364E-7*rho*(D_g*G_gOH+D_n*G_nOH)
POHm = 1.0364E-7*rho*(D_g*G_gOHm+D_n*G_nOHm)
PH2 = 1.0364E-7*rho*(D_g*G_gH2+D_n*G_nH2)
PH2O2 = 1.0364E-7*rho*(D_g*G_gH2O2+D_n*G_nH2O2)
PHO2 = 1.0364E-7*rho*(D_g*G_gHO2+D_n*G_nHO2)
PHO2m = 1.0364E-7*rho*(D_g*G_gHO2m+D_n*G_nHO2m)
POm = 1.0364E-7*rho*(D_g*G_gOm+D_n*G_nOm)
PO2 = 1.0364E-7*rho*(D_g*G_gO2+D_n*G_nO2)
PO2m = 1.0364E-7*rho*(D_g*G_gO2m+D_n*G_nO2m)
PO3m = 1.0364E-7*rho*(D_g*G_gO3m+D_n*G_nO3m)
PH2O = -(POH + POHm + 2*PH2O2 + 2*PHO2 + 2*PHO2m + POm + 2*PO2 + 2*PO2m + 3*PO3m)
PH = -2*PH2O-(PHp + POH + POHm + 2*PH2 + 2*PH2O2 + PHO2 + PHO2m)

H2O(0) = 4.110e+01
em(0) = 5.567e-09
Hp(0) = 2.746e-07
H(0) = 4.739e-09
OH(0) = 2.328e-11
OHm(0) = 8.635e-06
H2(0) = 1.480e-03
H2O2(0) = 4.903e-10
HO2(0) = 3.092e-14
HO2m(0) = 3.269e-12
Om(0) = 1.532e-13
O2(0) = 9.391e-16
O2m(0) = 2.200e-13
O3m(0) = 2.357e-25
rho = 0.73980
D_g = 0.0
D_n = 0.0
t(0) = 1e-07
t(f) = 25.000000000000007
